Cusps, self-organization, and absorbing states 
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Elastic interfaces embedded in (quenched) random media exhibit meta-stability and stick-slip 
dynamics. These non-trivial dynamical features have been shown to be associated with cusp singu- 
larities of the coarse-grained disorder correlator. Here we show that annealed systems with many 
absorbing states and a conservation law but no quenched disorder exhibit identical cusps. On the 
other hand, similar non-conserved systems in the directed percolation class, are also shown to ex- 
hibit cusps, but of a different type. These results are obtained both by a recent method to explicitly 
measure disorder correlators and by defining an alternative new protocol, inspired by self-organized 
criticality, which opens the door to easily accessible experimental realizations. 
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Elastic objects/interfaces evolving in disordered me- 
dia constitute a simple and unified way to describe phe- 
nomena as different as wetting, dislocations, cracks, do- 
main walls in magnets, or charge density waves In 
all these examples, two conflicting mechanisms compete: 
the tendency to minimize the elastic energy, leading to 
flat interfaces, and the propensity to accommodate to lo- 
cal minima of a random pinning potential, giving rise to 
rough profiles. In this way, many different meta-stable 
states exist, opening the door to the rich phenomenol- 
ogy characteristic of disordered systems: slow dynamics, 
creeping, stick-slip motion, avalanches, etc 0- At zero 
temperature and in the presence of a pulling force, a crit- 
ical point separates two different phases: i) a pinned one, 
in which disorder dominates and the interface remains 
trapped at some local minimum of the random potential, 
and ii) a depinned phase in which the force overcomes 
the pinning potential and the interface moves forward. A 
prominent example characterizing the dynamics of these 
systems is the Quenched Edwards- Wilkinson equation: 



dth{x, t) = iyV^h{x, t) + F + ri{x, /i(x, t)) 



(1) 



where v is a constant, h{x,t) the interface height, F 
an external force, and r]{x,h{x,t)) a quenched (random 
field) noise. The noise correlator is {r]{x,0)ri{x' ,u)) = 
Ao(u)5(x — x'), where Ao(u) is some fast-decaying func- 
tion. At a critical force, Fc, Eq.([T]) shows a depinning 
transition, representative of a broad universality class. 
Analytical understanding of Eq.(IT]) comes from various 
fronts, including the Junctional renormalization group 
(FRG) I, i, in which the disorder correlator, Ao(u), 
renormalizes as a function. This approach led, some years 
ago, to a successful computation of critical exponents in 
an e{— A — d) expansion up to one- loop order 3]. The 
emerging FRG fixed-point function for the disorder cor- 
relator, A*(m), presents the peculiarity of having a cusp 
singularity at the origin, i.e. A* (0"^) = —A* (0^) 7^ 0. 
Very interestingly, such "cusps" , rather than been a cu- 
riosity, have been argued to be essential to capture the 



physics of disordered systems 0, H, [E H- Physically, 
the underlying idea is that the effective energy landscape 
describing a random media in which an interface ad- 
vances consists of a series of parabolic wells (representing 
different metastable states) matching at singular points 
where the first derivative (i.e. the force) is discontinuous 
[3, H, 0I- In this way, the total random force experi- 
enced by the interface as it advances has necessarily a 
sawtooth profile, with linear increases followed by abrupt 
falls (see Fig.[T]). These jumps reflect the change from one 
metastable state to another; they dominate the statistics 
of the correlation functions and generate a cusp singular- 
ity at the origin of the renormalized disorder correlator. 
Le Doussal, Wiese, and collaborators have recently ex- 
tended the FRG calculation up to 2-loops 4] and devel- 
oped a strategy to measure disorder correlators, allowing 
to test the FRG predictions in computer simulations as 
follows [1,0. In Eq.((l|), F is replaced by a confining 
force Fldv/ = m^(w — h(x, t)), derivative of a parabolic 
potential centered at a, h = w, where w is a constant. 
For any w, a stable sample-dependent interface configu- 
ration, with average height h, is found. Then, by slowly 
increasing w, Fldw grows until, eventually, the interface 
overcomes a barrier and falls into a new meta-stable state 
with larger h (givingrise to sawtooth profiles). The main 
breakthrough in [1, Q is to prove that, for a size L, the 
cumulants oi w — h can be written as: 

{w~h{w)) ^ F,im)/m^, (2) 
{{w ~ h{w)){w' - Hw'))), = AUw - w')/im^L^), 
{[{w - hiw)){w' - Hw'))?)c = Sr^iw - w')/{m^L^''), 

and that, taking the limit m ^ 0, Fc{0'^) converges to the 
critical force Fc and Ao+ (resp. 50+) coincides with (is a 
function of) the FRG fixed point, A*{u) [1,0]. Numerical 
measurements of the cumulants in Eq.([5]) agree nicely 
with the 2-loop FRG predictions 

In this paper, we show that cusps appear also in the 
apparently very different realm of systems with annealed 
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disorder; systems with many absorbing states (AS), a 
conservation law, and no quenched disorder exhibit cusps 
identical to those of Eq. {Ip- Instead, similar models with 
AS in the directed percolation class, as well as the Bak- 
Tang-Wiesenfeld sandpile model have also cusps, but of 
different types. Moreover, we introduce an alternative 
strategy to measure cusps using self-organized critical- 
ity (SOC), i.e. by alternating slow-driving and boundary 
dissipation. This provides a practical and easy strategy 
to observe cusps in numerics and in experiments. 

Systems with absorbing states have annealed noise and 
a dynamics which leads to frozen/ absorbing microscopic 
configurations. Well known examples are directed perco- 
lation (DP), the contact process, or the Domany-Kinzel 
automaton ^8] . All these exhibit a transition from an ab- 
sorbing to an active phase in the very robust DP uni- 
versality class, represented by the Langevin equation: 
9tp(x, t) = ap — bp^ + V^p + ay^ri{x, t) where a and b are 
parameters, p(x, t) the activity field, and rj a zero-mean 
Gaussian white noise. The square-root noise ensures the 
AS condition: fluctuations cease at p = 0. Scaling fea- 
tures different from DP emerge only in the presence of ex- 
tra symmetries or conservation laws fs', 9| . Of particular 
interest are systems with many AS, in which absorbing 
configurations have a non-trivial structure, represented 
by a background field encoding the likeliness of absorbing 
configurations to propagate activity when perturbed, and 
allowing for metastability to appear. Two main classes of 
such systems are: 

(i) The Conserved-DP (C-DP) class (which captures 
the gist of stochastic sandpiles as the Manna or the Oslo 
one [id p is represented by a Langevin equation similar to 
that of DP but with an extra conservation law 0, [111, [13] : 

at/9(x,i) = ap - 6/9^ + 7p£:(x, i) + V V + 
dtE(ji,t) = DV^p. (3) 

a, b, D and 7 are parameters and i?(x, t) is the (conserved 
and non-diffusive) background field. Notwithstanding its 
similarity with the well-behaved DP equation, Eq. ^ has 
resisted all renormalization attempts and sound analyti- 
cal predictions are not yet available (IS] . 

(ii) The directed percolation class with many AS 
Defined by models as the pair contact process [3], this 
class has no extra symmetry /conservation-law with re- 
spect to DP. Its corresponding Langevin equation is: 

dtp{x,t) = ap-bp^ +-fp^{x,t) + V^p + ay/prj{x,t) 
at*(x,i) = ap-P^p. (4) 

where a, 6, 7, a and (3 are parameters. Despite the non- 
trivial absorbing phase, characterized by the background 
field *(x,t), Eq.lHD exhibits DP (bulk) criticality [U]. 

Remarkably, it has been conjectured that Eq.(IT]) and 
Eq. ^ are equivalent descriptions of the same underlying 
physics [ll|, . In spite of the absence of rigorous proof, 
there is strong theoretical and numerical evidence back- 
ing this conjecture [T^]. Indeed, there are heuristic ar- 
guments which, starting from Eq. ([3]) and defining a "vir- 
tual interface" H{x,t) = dt'p{x,t') which encodes the 
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FIG. 1: (Color online). Up (Down): Steady state time series 
of the spatially averaged force (background field) for a single 
run of Eq.ll} (Eq.Q) using the protocol by Le-Doussal-Wiese 
(using self-organized criticality) ; the average values have been 
subtracted in both cases. Linear increases are followed by 
abrupt falls: the fingerprint of cusps. 



past activity history, lead to Eq.([T|) as an effective equa- 
tion for H{x,t). Note that any past noise trajectory is 
a frozen variable, allowing to map present-time annealed 
noise into quenched disorder. Moreover, the background 
field E{x,t) in Eq.Q can be identified with the force 
F + i^V^h{x,t) in Eq.((T]); both encode the propensity to 
propagate activity/motion at each site p^ . 

In what follows, we a) reproduce the results in Q for 
Eq.([T]) and introduce an alternative protocol to measure 
them, b) look for cusps in Eq.Q and compare them 
with those in quenched systems to verify if, indeed, they 
both represent the same physics, and c) explore whether 
different types of cusps exist for other AS systems. 

a) We directly integrate a discretized version of Eq. ([T]) 
in one dimension, with a force m^{w — h(x,t)) and pe- 
riodic boundary conditions 0. r]{'x.,h) takes quenched 
random values at discrete equispaced values of h, and 
is linearly interpolated in between. To simplify the nu- 
merics we use an overall Heaviside step function in the 
r.h.s. of Eq.lJT]) forbidding the interface to move back- 
wards [16] . For each w the dynamics eventually reaches 
a pinned configuration (non-positive r.h.s. everywhere). 
By increasing quasi-statically the value of w and using 
Eq.([2]), we determine Fc, A(u) and S{u). It is conve- 
nient (see 5, 6]) to use normalized functions defined by 
Y{u/u^) = A(u)/A(0) with fixed by J^Y{z)dz = 1 
and Q(A(ii;)/A(0)) = J^' Siw')dw' / Siw')dw' . The 
universal functions Y{z) and Q{y) are determined for up 
to L = 2^^ and = 0.001. Results are indistinguishable 
from those in Q for both functions (see Fig. [2). 

We now define an alternative protocol to measure 
cusps. The basic idea is to let the system self-organize 
to its critical state by iterating local slow driving and 
boundary dissipation as done in self-organized criticality 
[13, [IB]- Starting from a pinned interface, slow driv- 
ing is implemented by increasing the force F (which be- 
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FIG. 2: (Color online). Upper(Lower): Y{z) {Q{y)) for nine 
different one-dimensional cases: i) Ec|Jl| for random field 
and ii) random bond disorder (from and iii) from our 
own simulations, iv) self-organized Eq.([T]|, v) C-DP, Eq.((3}, 
vi) self-organized Eq.Q, vii) self-organized Oslo model and, 
finally, viii) a critical conserved reaction diffusion model. All 
these curves collapse into a unique one, different from the one 
obtained for ix) DP with many AS, Eq.Q, following the self- 
organized procedure. Differences are larger for Q{y) than for 
Y{z). Inset: function S{u), from which Q{y) is determined, 
for Eq.©. 



comes a x-dependent function) at a randomly selected 
site by a small amount. This is repeated until the in- 
terface gets depinned at some point. Then (and only 
then) the dynamics, Eq.([T]), proceeds. Open boundaries 
h{0) — h{L -I- 1) = allow for "dissipation" or force- 
reduction (i.e. owing to the open-boundaries, the inter- 
face develops, on average, a parabolic profile with neg- 
ative average elastic force) inducing the system to get 
pinned again, and so forth. This driving/dissipation cycle 
is iterated until a steady state is reached. Then correla- 
tions between the spatially-averaged (pulling plus elastic) 
force F + V^/i, at different driving steps, are measured. 
Note that this averaged force defines a stochastic (gen- 
eralized Ornstein-Uhlenbeck) process confined inside a 
L-dependent potential. As the force balances at each 
pinned state the corresponding random pinning force, 
its correlations provide direct access to the quenched- 
disorder correlations. We have verified that the aver- 
age pulling force is identical to the critical force Fc (see 
Eq.(I21)) and measured the cumulant ratios; both Y{z) 
and Q{y) turn out to be indistinguishable from their 
counterparts above (see Fig. ^ . This coincidence follows 
from the fact that we are just considering two alternative 
"ensembles" to study Eq.([T|) at criticality. In the first. 



fixed velocity ensemble^ one considers closed boundaries 
and lets the velocity go to 0+ (m^ 0). Instead, in the 
second one, we let the system self-organize by combining 
slow driving and boundary dissipation at infinitely sep- 
arated timescales fisl. These two paths lead to the 
same critical properties [l9l |. 

b) We now adapt these two protocols to study sys- 
tems with many AS. We start analyzing Eq.Q in one- 
dimensional lattices (size up to L = 2^^) using the re- 
cently proposed efficient integration scheme for Langevin 
eqs. with square-root noise [20|. The system is initialized 
with arbitrary initial conditions and then it self-organizes 
to its critical point by means of slow driving and bound- 
ary dissipation (second protocol) as follows: When 
the system is at an AS, the background and the activity 
field are increased, at a randomly chosen site, by a small 
amount {SE — Sp = 0.1); then, the dynamics proceeds 
according to Eq.Q. Open boundaries allow to decrease 
the averaged background field and hence the linear term 
in the activity equation, until eventually a new AS is 
reached. Then, a new driving step is performed and so 
on, until a steady state is reached. Fig. [T] shows the evo- 
lution of E (minus its average value) as a function of the 
number of driving steps in such a state; linear increases 
are followed by abrupt decays corresponding to dissipa- 
tive avalanches. From such time series, we measure the 
three cumulants of the spatially averaged background 
field (equivalent to the force above) and, from them, the 
critical point, Y{z) and Q{y)- Results for L = 2^°, sum- 
marized in Fig. [21 show an excellent agreement with their 
counterparts obtained for Eq.([T|). 

To obtain results for Eq.([3|) using the strategy in 
one needs to introduce a slowly moving parabolic poten- 
tial {m^{w - i?(x,<)) with -ff(x,i) dt' p{x,t')) for 
the background field and periodic boundaries. This is 
achieved by including a term ^m? p['x., t) in its Langevin 
Eq.(l3]) and increasing i? by a small constant amount, 
rri^Sw, after each avalanche. When integrated in time, 
these two contributions give a slowly moving parabolic 
potential for E. Additionally, after each avalanche, ac- 
tivity at some random point is created to avoid remaining 
trapped in the AS. This causes a cascade of rearrange- 
ments in the activity and background fields which even- 
tually stops owing to the dissipation in E induced by 
the negative forcing term {—m^p). From the correlations 
of average background field values in the AS at different 
steps, we obtain results identical to those reported before 
for both Y{z) and Q{y) (see Fig. [2]). 

The (impressive) overlap of whole functions makes a 
much stronger case for the conjectured equivalence be- 
tween Eq.([T]) and Eq.Q (and also between the two pro- 
tocols) than any previous numerical agreement of critical 
exponents. These conclusions have also been extended to 
bidimensional systems (not shown). 

Let us now study different stochastic sandpile mod- 
els [lo| as well as a reaction diffusion system with con- 
servation 0, all of them argued to be in the C-DP 
class [31 but for which, owing to the discrete nature of 
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particles/grains, the continuous statistical tilt symmetry 
(which leaves the physics of Eq.(IT]) invariant under the 
continuous h{x) h{x) + 5h{x) transformation Q) is 
broken. Using the self-organized protocol, we produce 
the cusps reported in Fig. [2] The curves overlap al- 
most perfectly with the previously obtained ones, im- 
plying that the continuous statistical tilt symmetry is 
asymptotically restored and confirming the universal be- 
havior of stochastic sandpiles and their associated cusps 
[21I I . Instead, for the deterministic Bak-Tang-Wiesenfeld 
sandpile [13], known to exhibit different critical behav- 
ior, we have measured different correlators, with another 
type of cusps (not shown). 

c) Finally, we apply the previous methods to the DP 
(with many AS) class, Eq.([4]), for which an effective in- 
terfacial description can also be constructed [1^ . Results 
are summarized again in Fig. [2l Y{z) shows a cusp at 
the origin which differs slightly from the one above; the 
discrepancy in Q{y) is much more pronounced, confirm- 
ing that DP exhibits a different type of scaling, and that 
cusps are a common trait of systems with many AS. This 
is a remarkable, so far unveiled, property of DP systems. 

In summary, we have shown the presence of identical 
cusps in systems with quenched disorder, systems with 
AS and a conservation law, and self-organized stochas- 
tic sandpiles. In contrast, different cusps are obtained 
for systems in the directed percolation class as well 



as for the deterministic Bak-Tang-Wiesenfeld sandpile. 
This enriches our present view of universality in stick- 
slip/avalanching systems and confirms the ubiquitous 
presence of disorder correlators with cusps in avalanching 
systems, with or without quenched disorder. 

Our new protocol to determine cusps opens a straight- 
forward path to measure them in experiments. Note that 
for a sandpile it suffices to measure its weight after every 
avalanche to determine cusp- correlators! (see [1^ for a 
student-laboratory set-up). Experimental studies of the 
mass of granular piles have actually noticed the presence 
of sawtooth fluctuations (see Fig. 1 of [13]) and, for exper- 
imental ricepiles, their correlations have been looked-at 
[2^, but without determining the relevant ones. Simi- 
larly, cusps could be easily measured in any system, as 
superconductors [27j . in which self-organized criticality 
has been observed in the laboratory. 

Challenging questions remain unanswered: Could the 
FRG cusps be derived from a renormalization group so- 
lution of Eq.(l3])? Could the DP cusps be determined an- 
alytically from Eq.Q? Are there further systems that 
exhibit cusps? Answering these questions and study- 
ing experimental realizations, would greatly enhance our 
understanding of stick-slip/avalanching dynamics, and 
would provide a rich cross-fertilization between interfaces 
in random media and systems with AS. 
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